Existence threshold for the ac-driven damped nonlinear Schrodinger solitons 
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It has been known for some time that solitons of the externally driven, damped nonlinear 
Schrodinger equation can only exist if the driver's strength, h, exceeds approximately (2/77)7, where 
7 is the dissipation coefficient. Although this perturbative result was expected to be correct only 
to the leading order in 7, recent studies have demonstrated that the formula htbx = (2/77)7 gives a 
remarkably accurate description of the soliton's existence threshold prompting suggestions that it 
is, in fact, exact. In this note we evaluate the next order in the expansion of /ithr (7) showing that 
the actual reason for this phenomenon is simply that the next-order coefficient is anomalously small: 
hthr = (2/77)7 + O.OO27 . Our approach is based on a singular perturbation expansion of the soliton 
near the turning point; it allows to evaluate /ithr (7) to all orders in 7 and can be easily reformulated 
for other perturbed soliton equations. 

PACS number(s): 03.40.Kf, 42.65.Tg, 42. 81. Dp 



I. INTRODUCTION 

The externally driven, damped nonlinear Schrodinger 
(NLS) equation, 



iStft + * M + 2|*| 2 * = -27* - he 



'in 



(1) 



arises in a variety of fields including plasma and con- 
densed matter physics, nonlinear optics and supercon- 
ducting electronics. In some of these applications (e.g. in 
the study of the optical soliton propagation in a diffrac- 
tive or dispersive ring cavity in the presence of an input 
forcing beam JjJ ; in the description of easy-axis ferromag- 
nets in an external rotating nragnetic field perpendicular 
to the easy axis JH]; in the theory of rf-driven waves in 
plasma ||) Eq.(|l]) has a direct interpretation. In oth- 
ers — like for instance in charge-density- wave conductors 
with external electric field [ 1 ; shear flows in nematic liq- 
uid crystals H; ac-driven long Josephson junctions H, 
and periodically forced Frenkel-Kontorova chains fij] 
it occurs as an amplitude equation for small and slowly 
changing solutions of the externally driven, damped sine- 
Gordon equation: 



qtt + Mt 



sing = Tcos(wi). 



Without loss of generality f2 in Eq.p) can be normal- 
ized to unity [p|-|lO||; hence, the driver's strength h and 
dissipation coefficient 7 are the only two essential con- 
trol parameters. Given some h and 7, a fundamental 
question is what nonlinear attractors will arise at this 
point of the (7, /i)-plane. In their pioneering paper |ll| | 
Kaup and Newell considered Eq.(fj) on the infinite line 
under the vanishing boundary conditions at infinity. By 
means of the Inverse Scattering-based perturbation the- 
ory, these authors have demonstrated that for small h 



and 7 Eq.(p]) exhibits two soliton solutions phase- locked 
to the frequency of the driver. As h is decreased for 
the fixed 7, the two solitons approach each other and 
eventually merge in a turning point for h = (2/77)7 PI - 
Consequently, this value plays the role of a threshold; no 
solitons exist below h — (2/77)7. Later the same exis- 
tence threshold was reobtained by Terrones, McLaugh- 
lin, Overman and Pearlstein H in a regular perturbation 
construction of solutions to (|l]) in powers of h and 7 (see 
also @). 

In ref. J10[ equation (|lj) was studied, numerically, in 
the full range of h and 7. It was found that the two soli- 
ton solutions persist for 7 up to approximately 0.7. For 
each 7 < 0.7 there is a turning point at some h — hthr 
at which one branch of solitons turns into another, and 
which plays the role of the lower boundary of the exis- 
tence region |lj] . Amazingly, Kaup and Newell's approx- 
imate relation /i t hr = (2/77)7 was found to remain valid 
even for not very small 7. For example, for 7 = 0.48 the 
ratio /ithr/7 is different from 2/77 by only one part in a 
thousand M). 

A completely different approach was put forward b 
Kollmann, Capel and Bountis |l4J who regarded Eq. 
as the continuous limit of a discrete NLS equation which 
they studied by means of the fixed point analysis and 
the Mclnikov-function method. In particular, the lower 
boundary was obtained from the tangential intersection 
of the invariant manifolds of a hyperbolic fixed point. A 
renrarkable accuracy of Kaup and Newell's linear law de- 
tected in ]Hj as well as conclusions of their own Melnikov- 
function analysis prompted the authors of Q to suggest 
that the relation /i t hr = (2/77)7 can be exact, at least for 
sufficiently small 7. 

The aim of the present note is to demonstrate that this 
relation is, in fact, not exact, and the actual reason why 
it appears to be so accurate for small 7 is simply because 



the coefficient of the next term in the expansion of hthr(j) 
in powers of 7 is anomalously small. We do this by re- 
constructing the two solitons in the vicinity of the lower 
boundary of their existence domain by means of a singu- 
lar (rather than regular) perturbation expansion. This 
novel expansion constitutes the main technical achieve- 
ment of our work; its scope of applicability is obviously 
much wider than just the damped-driven NLS equation 

The reason why the regular expansion is not adequate 
near the solitons' existence boundary is quite simple. 
The point is that the regular expansion is based on the 
assumption that solutions can be continued along rays 
h = f)7 (f) = const.) But since the boundary is not ex- 
actly a straight line h = (2/77)7 ( as wn l be shown below, 
it slowly recedes upwards from this straight line), the ray 
h = f)7 with t) slightly above 2/V will hit the boundary at 
some small 7 = 70. Consequently, the regular expansion 
fails to continue the solution beyond 70- In the singular 
expansion, on the other hand, the solution is continued 
along a curve whose shape is calculated simultaneously 
with finding perturbative corrections to the solution it- 
self. In this way the existence boundary can be found 
to any desirable accuracy. (In this paper we restrict our- 
selves to terms ~ 7 3 .) This idea should remain applicable 
to other perturbed soliton-bearing equations. 

The outline of this paper is as follows. We start by dis- 
cussing the regular asymptotic expansion as h and 7 — > 
(section O) . The procedure is similar to the one in M ; the 
only difference is that since we now deal with solutions 
decaying at infinities (^> x — ► 0) rather than periodic as 
in ||], we will be able to find perturbative corrections in 
closed form. In section III we explain why the perturba- 



tion series for ^f breaks down as h approaches the turning 
point, and replace it by a singular expansion. This allows 
us to find the next terms in the expansion of /ithr(7)- The 
resulting asymptotic formula is compared then with the 
threshold /tthrW obtained in a high-precision numerical 
analysis of Eq.(g) for several values of 7. Next, after we 
have achieved an understanding of why the regular ex- 
pansion fails and how it can be cured near the existence 
threshold, a natural step is to try to develop a unified ex- 
pansion which would be equally applicable near and far 
from the threshold. This is done in section pM Finally, 
some concluding remarks are made in sectionjV] followed 
by a brief summary of our results. 



II. REGULAR PERTURBATION EXPANSION 

By making a substitution ^f(x, t) — ip(x,t)e lt Eq.(|l]) 
can be reduced to an autonomous equation 



iipt + ipxx + 2|i/'| 2 V' -ip = -ijijj - h. 



(2) 



We will be interested in time-independent solutions of 
Eq.(0); these satisfy the stationary equation 



with the boundary conditions 

ipx(%) ^0 as \x\ — ► 00. 

We start with developing a regular perturbation ex- 
pansion away from the turning point. As the authors of 
M , we assume that we are approaching the origin on the 
(7, ft)-plane along a straight line h = F)7 where f) is a 
proportionality coefficient. Letting 



ip=(u + iv)e- ia , 



(4a) 



where a is some constant phase that can be conveniently 
chosen at a later stage, we expand 



u = Uq + 7«i + ..., v = v + 7V1 + 



(4b) 



and substitute into Eq. (||) . The coefficient of 7 gives the 
unperturbed stationary NLS equation with a well-known 
soliton solution 






cos 9 



sechx. 



Here 9 is a free parameter. Next, at the order 0(7 1 ) one 
gets 



H 



»1 



f) cos a — Vq 



V\ j \ \) sin a + Uq 
where the Hermitean operator 



(5) 



"■- !-l^+iV-2fo?,+ 3M ° ..,rT..2 s ^ 



dx 2 



Vq + 3mq 2u v 
2u v u^ + 3v 2 



and I is the 2x2 identity matrix. In order for the equa- 
tion (||) to be solvable in the class of bounded functions, 
its right-hand side needs to be orthogonal to the vector 
(vo, —uo) T , the eigenfunction of the operator Hq associ- 
ated with the zero eigenvalue. (This zero eigenvalue re- 
sults from the U(l) phase- invariance of the unperturbed 
NLS equation.) The orthogonality gives a relation be- 
tween a and 9, 



7rhsin(0 — a) = 2, 



(7) 



tp xx + 2\ip\ 2 ip - ip = -iji() - h 



(3) 



implying that only one of the two parameters (say, 6) 
can be chosen freely. It does not matter what exactly 
we choose for 9; the net phase of the leading-order 
approximation is (9 — a) and this is fixed by Eq.(R). 
The meaning of this relation is straightforward. For 
h = 7 = 0, the NLS equation has a family of soliton 
solutions, ip — e l ( ~ a * > sechx, with (6 — a) arbitrary How- 
ever, if we want to continue the solution along the line 
h = f)7, the unperturbed solution that we need to start 
with has the phase given by Eq.(R). 

It is convenient to take 9 — 7r/2; this makes the lin- 
ear operator (ra) diagonal. (The other diagonal choice 
9 = is also acceptable but somewhat less convenient 



in the present context.) The constant phase a is then 
determined by 



2 1 
cos a = — . 

7T \) 



(8) 



In fact, there are two values of a defined by this equa- 
tion, one positive and one negative. The positive a = a + 
corresponds to the soliton ijj^ + > and the negative a = a_ 
defines the soliton ijj(~\ Since the left-hand side cannot 
exceed 1, the right-hand side gives the well-known for- 
mula for the lower boundary of the domain of existence 
of the two solitons: t) > f) t hr = 2/tt (H^JH. (In the 
next section we will obtain a more precise formula for 
this threshold.) 

Now for 9 = 7r/2 the equations (ph become 



LqU\ (x) — f) cos a — vq (x) ; 
L\Vi(x) — f) sin a, 



(9) 
(10) 



where Vo(x) = sechx and L$ and L\ are the standard 
Schrodinger operators with familiar spectral properties: 



L = -d 2 /dx 2 + 1 - 2sech 2 a;; 
L x = -d 2 /dx 2 + 1 - 6sech 2 x. 



(11) 
(12) 



The operator L\ is invertible on even functions; in par- 
ticular, 



L 1 "^echa; = -(xtanhx— l)sechx, 



L, sech x = — sechx 
4 



and 



(13) 



(14) 



(15) 



L^ 1 1= l-2sech 2 :r. 
Hence Eq.(|10|) is readily solved: 

Vi = f)sina(l — 2sech 2 ir). (16) 

The condition (g) being in place, Eq.(g) is solved as well: 



u\ = U\(x) + Asechx, 



(17a) 



where 



2 1 1 

U\{x) = — I — tanhxsinha; -\ — x {j(x) sechx — 

7T 2 7T 

— (xsechx + sinhx) arcsin(tanhx) — 1}, (17b) 

j(x) = / £sech£ci£, 
Jo 

and A is an arbitrary constant which is to be fixed at 
higher orders of the expansion. Hence we proceed to 
0( 7 2 ) to find 



L u 2 = (4:V ui - l)ui, 
Lxv% = Mi + 2v a (u\ + iv\). 



(18) 
(19) 



Equation ([FJ) is solvable if its right-h and side is orthogo- 
nal to sechx. Substituting from (|l6|)-(17a), this condition 
fixes the constant A: 



AW = 1 



4 f n , 

— / Wi(a;)sech x(l — 2sech x)dx, 

71" J 



(20) 



where we have used Eqs.( |l4| - |l5| ). (Here we have written 
A^ for A so as to emphasize that this is now a fixed num- 
ber; this number will reappear in the singular expansion 
below.) Eq.(p8[) is now solved in the form 



u 2 = Uzix) + Bsechx. 



(21) 



The constant B is to be fixed at the 7 3 -level, where we 
obtain the equation 

L Q u 3 = 2{ui(u 2 + u\) + 2v (u 2 vi + uiv 2 )} - v 2 . (22) 

The solvability condition for eq.(|2^) gives us B: 

B = (7r()sina) _1 x 

vq{2ui{u\ + v 2 ) + Avq(uiv 2 +U 2 vi) — v 2 }dx. (23) 



So far our treatment followed the lines of Terrones et al 
P]; the only difference is th at ou r vq,Ui,Vi, ... are given 
by explicit formulas. Using (17b) in ( poj ) and integrating 
numerically, we identify the constant A*- ' which com- 
pletes the determination of the first-order corrections: 
A@) = -2.4378 x 10" 1 . 

Let us now send h — ► (2/71)7. The formula (O) for 
U\{x) is not affected and the expression (|20| ) for A re- 
mains valid as well. Therefore, the solvability of Eq.(fi~|) 
is ensured and u 2 can be written in the form (|2lj). The 
constant B is expected to be identifiable from Eq.([23|). 
However, for F) — » 2/7r we have sin a — > and so this 
formula gives B = 00 unless 



vq(2u\ + ivoUiv 2 — v 2 )dx = 0. 



(24) 



(Here we have used that v% — > as sin a — » 0.) In gen- 
eral the condition ( |24| ) is not in place, and therefore the 
regular expansion blows up. 



III. SINGULAR PERTURBATION EXPANSION 
AT THE TURNING POINT 

The reason for the breaking down of the expansion is 
that it was implicitly assumed in Eq.(4b) that vi — 0(1) 
whereas in the actual fact, in the limit t) — > 2/7T we have 
vi — ► 0. Let us now explicitly take this fact into account 
by writing 



u = jm + 7 u 2 . 



Vq + 7 v 2 + 



(25a) 



where vq = secha;. We also expand \): 



f) = f)Q + f)l7 + ()27 



fa 



(25b) 



(Thus we have fixed 8 = tt/2 and a — straight away.) 
Substituting into (f|), the first order in 7 yields eq.(Q) 
where one should only replace t| cos a — > fa. Its solu- 
tion is given by the same eq.(|17]) as before, with A an 
undetermined constant. At the order 7 2 we obtain 

L u 2 = fa, 
and hence \)i = and u 2 = B secha;. The equation for v 2 



L\v 2 = u\ + 2vqu 1 



(26) 



[cf. Eq.(J19|)1; this is always solvable. Finally, the 7 3 -level 
yields 

L u 3 = 2(ul + 2v Q u 1 v 2 ) - v 2 + fa 
[cf . (122)] , whose solvability condition is given by 



vq (2u\ + 4v uiv 2 — v% + fa) dx = 0. 



(27) 



It turns out that it is on ly this equation ( |27J ) that fixes 
t he c onstant A in E q.(|l7a| ). Indeed, substituting u\ from 
( |!7a| ) and v 2 from (|26j), Eq.(p7|) reduces to a quadratic 
equation 

A 2 - 2PA + Q - nfa = 0, (28) 

where after some algebra the coefficients are found to be 

P = -2t ) l + ^-J{fa-U 1 (x)}dx (29) 

and 

Q = 1{U\ - Wi(l + 2Wisechx) x 
x [f) (l — 2sech 2 a;) + secha; (1 — x tanha;)]}da;. (30) 

In the derivation of (E9-p0h we used Eq. (M) and the iden- 
tity 



4L^ 1 (v 2 U 1 )=L^ 1 (t)o-v )-U 1 ; 



(31) 



this is a straightforward consequence of Eqs.(H),d9|) and 
the fact that the Schrodinger operators ( 11 )-([12[) differ 
by 4sech 2 a;: 

L = L 1 +Avl{x). 

Since there is a cubic term in u\ in Eq.(|27j), one could 
expect the resulting equation for A to be cubic; however 
the coefficient in front of A 3 is easily shown to vanish. 
Another observation is that the coefficient P coincides 
with the constant A^ [Eq. (pOJ)] obtained in the regular 



expansion. To see that, one only needs to use the identity 
(pl| ) once again. 

The roots of Eq. (p8|) are given by 



A (±) = A (o) ± n U _ f, 



,(o) 



where 



I) 



(o) _ 



-XQ-P 2 )- 



(32) 



(33) 



Hence, if \] 2 > fa 2 , Eq.(H) has two solutions i/j^ which 
are only different in the coefficients A^ 



If fa < t)i°\ 



(0) 



there are no solutions at all. The value fa = t) 2 * s there- 
fore the turning point. Doing numerically the integral in 
@ we find Q = 6.4665 x 10" 2 . Recalling that P coin- 
cides with Eq.©, P = A(°) = -2.4378 x 10~\ Eq.(||) 

gives fa — 1.667 x 1CP 3 . Finally, the coefficient A cor- 
responding to the turning point coincides with the "far 
from the turning point" value, Eq.(pCJ): A = A^°\ 

It is worth noting here that if fa = 0, Eq.(g7|) i s f° r_ 
mally coincident with Eq. (EJ) . This does not mean, how- 
ever, that soliton solutions exist for t)2 = 0, and that 
these solutions can be found by regular expansions (0) . 
The difference is that in Eq.(|24j) the function u\{x) has 
the coefficient A which has already been fixed by Eq. (EOh , 
whereas Eq.(27) is an equation for unknown A. 

Next, how close to hthx does the regular expansion stop 
working and has to be replaced by the singular one? For 
f) = 2 /it + f)27 2 Eq.(|) produces 

a± = ±v/^ 7 + 0(7 3 ) (34) 

and so the regular expansion for ip^' reads 

tp^' = i secha; + jUi (x) + 
+7 (A {0) ± y/^fa] secha; + 0( 7 2 ) (35) 



with U\ as in (17b), whereas the correct, singular expan- 
sion is 

^(±) = i secha; + -fUi(x) + 7v4 (±) secha; + 0(7 2 ), (36) 

where A^ are given by Eq.(|32|). Comparing ( p5|) to 
(pq), one concludes that the difference between the reg- 
ular and singular expansions is negligible provided l)2> 
f^ 0) w 2 x 10~ 3 . On the contrary, for f) 2 ~ 10~ 3 the 
difference cannot be ignored. 

We close this section by comparing our asymptotic ap- 
proximation for hthr (7) with results of the high-accuracy 
numerical solution of Eq. (|3|) . Here we employed a fourth- 
order Newtonian algorithm, with the stepsize Aa; = 0.005 
and residual 5 x 10 -11 . Eq.(j3|) was solved on a half- 
interval (0,40). We would decrease h in increments of 
Ah = 1 x 10~ 10 until iterations have ceased to converge; 
after that we would repeat the computations on a con- 
densed grid and extended interval. As a result we were 



able to obtain the threshold h t \ 1Y accurate to ten digits 
after the decimal point. 

Fig.l shows the difference between the numerically 
found values of ft-thr and (a) Kaup and Newell 's linear 
law /ithr = (2/ 7r )7i an d (b) the refined asymptotic ex- 
pansion hthr = (2/7t) 7 + ^ 0) 7 3 with t) { 2 0) = 1.667 x 1(T 3 . 
The latter expansion is indeed seen to provide a more 
accurate approximation for the threshold. 
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FIG. 1. The comparison of the numerically ob- 

tained threshold values and their analytic approximations. 
Plotted is the difference (ftthr(n«mericoi) ~ h thl(analytic) ) /-y, 
where h thT{analytic) is given by (2/77)7 (dashed curve) and 
(2/tt)7 + 0.00177 3 (solid line). 



IV. A UNIFIED VIEWPOINT 

The singular expansion of the previous section was de- 
signed as a continuation and generalization of the regu- 
lar expansion and presented in the form that allows for 
a straightforward comparison with the latter. A natural 
step now is to try to develop a unified approach which 
would be valid both near and far from the turning point. 
Such a unified formalism could provide an additional in- 
sight into the structure of solutions and have some tech- 
nical advantages. 

The unification is achieved if one notices that the 
role of the coefficients A, B, ... is simply to renormalize 
the constant phase of the leading-order approximation 
■00 = isechx. Consequently, instead of adding homoge- 
neous solutions Asechx, Bsechx, ... at each order of 7 we 
can expand a: 



a + ai7 + a 2 7 + 



(37) 



(This observation belongs to D. E. Pclinovsky.) Multi- 
plying Eq.(0) through by e la : 



ip xx + 2|V>| 2 V' ~ip = -h^P - he 1 



and substituting Eqs.(25b), (37) and ( Bq ) with uq 
and vq — sechx, we get, at the order 7: 



(38) 
= 



2 
f)o = — 



1 



7T cos ao 



(39) 
(40) 



with U\ as in Eq.( |l7b|) . For any ao Eq.(|39|) gives the 
initial slope of the curve h(j) along which we want to 
continue our solution. The lowest h arises for ao = 0; 
hence in order to obtain the threshold one has to set 
()o = 2/ir. Next, the order 7 2 yields 

fji = fjotana |aH — / v (l - 4v a ui)L^ 1 ldx > (41) 



and 



v 2 = L x 1 (ui + 2v ul + t) ai). 



(42) 



Eq.(fiT]) relates f)i and a.\. One of these (say, f)i) can be 
chosen freely; then Eq.(pd|) fixes the other. This simply 
means that given a nonzero ao , the solution exists for an 
arbitrary \\\ (i.e. in an arbitrary neighbourhood of the 
ray h — f)o7-) The case ao = is exceptional; in this 
case the ray is very close to the threshold and so we have 
to set f)i = while ai remains undetermined. (This is 
the case where we would have to invoke the singular ex- 
pansion before.) The coefficient ai can only be fixed at 
the 7 3 -level, where we obtain 

- v {2[ui(ul + v\) + 2vo(u 1 v 2 + u 2 vi)] -v 2 }dx = 



< r r / ■ U 1 

= f]iai sin ao — 1J2 cos a o + 1)0 I ot.% sin ao + — cos ao 



Letting ao = this becomes exactly the quadratic equa- 
tion (Eq) (where we only need to replace A — > ol\) and 
hence we have reproduced our previous threshold. 

Note that in the new formalism the corrections in the 
phase of ipo are uncoupled from the rest of the expansion. 
The advantage of this is that the correction ai (= A in 
the previous notation) arises at the order j 2 and not at 
0(7) as before; the correction a 2 (previously known as 
B) arises at 0(7 3 ) and not at 0(7 2 ), and so forth. As 
a result, u\ does not contain the unknown ai (A); the 
correction v 2 depends on ai only linearly (not quadrat- 
ically as before), and so the resulting equation for ai is 
manifestly quadratic. (Recall that Eq.(E8h arises initially 
as a cubic equation and only then the coefficient of the 
term A 3 is calculated to be zero.) 



V. CONCLUDING REMARKS AND 
CONCLUSIONS 



1. In the undamped case (7 = 0) for any h e (0, ^/2/27) 
Eq.(f|) has two explicit solutions g: 



?i (±) W^Ji + 



2 sinh 2 (3 



1 ± cosh(„4x) cosh (3 J ' 



(43) 
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(all three pairs of solutions are 
practically indistinguishable) 



FIG. 2. Soliton transformation for small 7 (7 = 0.01). 
(a,b): 

For h far above the turning point h t hi = 6.3661997 x 10 3 
the imaginary parts are close to zero and the two solitons are 
well approximated by the two undamped solitons (W3J) (dashed 
lines). As h approaches the turning point, the real parts de- 
crease and converge while imaginary parts grow (dotted then 
solid lines.) (c,d): In the immediate vicinity of /i t h r the hump 
in the tp^ + ' profile rapidly transforms into a dip and the two 
solitons collapse into one (dashed then solid then dotted line.) 
All functions being even, we only show them for positive x. 



where 



^ = {2(1 + 2 cosh 2 /?)} V2 



is the asymptotic value of ip >{x) as |x| — > 00; the pa- 
rameter P is defined by inverting the relation 



h = 



\/2cosh 2 /3 
(l + 2cosh 2 /3) 3 / 2t 



and A is given by 



A — 2-0 oo sinh /3 



V2sinh/3 
v / l + 2cosh 2 /3' 



For the given value of h the tp^ branch merges 
with -i/A - ) at some 7thr [defined, approximately, by h = 

(2/7r)7thr + f)2 7thJ- If ^ is sm all then this 7thr is also 
small, so that the point of the merger is very close to the 
/i-axis and hence intuitively one could expect solutions 
at this point to be close to the undamped solitons (|43|). 
However, it is not quite obvious how this proximity can 
be reconciled with the fact that the two real solutions 
( p3| ) are rather far from each other. 

A related question concerns the shape of the ip( +s> and 
■0(~) solitons. In the undamped case [Eq.(43)] the real 
function -0( + ' has a hump and ip(~> has a dip. The nu- 
merical analysis shows that the hump respectively the 
dip persist in real parts of ■0 ( - + ' respectively ip(~~> soli- 
tons for 7 7^ |l(J. Again, it is intuitively not quite 
obvious how the hump can transform into dip as the two 
branches merge. Shall this transformation proceed via 
the flat state? 

The behaviour of solutions for small h and 7 (in partic- 
ular, near the turning point) can be clarified by invoking 
the asymptotic expansions (||) and (|2^). For the sake 
of illustration, we have also computed the two solitons 
numerically for a fixed small 7 (7 = 0.01) and varying 
h; results are shown in Fig. 2. In agreement with Eq.([|), 
for not very small h (h > 0.1) the imaginary parts of 
the solitons are seen to be almost zero [dashed lines in 
Fig. 2(b)] while the real parts change slowly with the vari- 
ation of h. For these h the pair of real solutions (H3) 
does indeed provide a good approximation for the corre- 
sponding -0W with 7 = 0.01. However, as h goes down, 
the real parts start changing (decreasing in magnitude) 



more vigorously while the imaginary parts are not small 
any longer; consequently, the approximation deteriorates. 
Near the turning point the real parts of both solitons be- 
come much smaller than their imaginary parts [see Fig. 
2(c,d)]. Nevertheless, for h not very close to the turning 
point (more specifically, for h > 6.3665 x 10~ 3 ), the real 
part of ip^ still has a hump and real part of i)A^^ still 
has a dip [Fig. 2(a)]. This justifies our usage of the no- 
tations ifj( + > and ip(~> for 7 7^ 0. Finally, in a very near 
vicinity of the existence threshold hthi = 6.3662 x 10~ 3 , 
the hump of the real part quickly transforms into the dip 
[Fig.2(c)]. _ 

This rapid change near the turning point can be eas- 
ily understood in terms of the singular expansion (25) 

When 1)2 is close to t} 2 
ing 



(0) 



we can define a small e by writ- 



7Tf)2 = 7Tt)2 + e- (44) 

The corresponding A's in Eq.( |l7a| ) are then given by 

A {±) = A (o) ± ^ (45) 

Recalling that to the leading order in 7 it is only this co- 
efficient A that determines the dependence of solutions 
on h, Eq.(p5|) gives the rate of change of their real parts: 



— = ±-J- S echx + 0( 1 2 . 
d\) 2 2 Ve 



(46) 



As e — > 0, the rate of change becomes infinitely large. 
Away from the neighbourhood of the turning point, in 
the region of the applicability of the regular expansion 
(0) , the rate of the transformation of the solitons ip± 1S 
given by 



da 



1 



hthr = -7+ (1.667 x 1(T 3 )7 3 + 

7T 



(47) 



(b.) For h away from the above threshold, more pre- 
cisely for h — (2/-7r)7 3> O.OO27 3 , the solitons are given 



where vi and u\ 
-@ with A<® = 



by the asymptotic expansion Eq.(0 
are given by explicit expressions (|l( 
-2.4378 x 10" 1 . 

(c.) For h close to the turning point, h = (2/ir)j+i)2j 3 
with t)2 ~ 10~ 3 the second, 7 2 -order of the regular ex- 
pansion (||) becomes greater than the first order, and the 
expansion breaks down. In this case the two solitons 
are given by the singular expansion (|36|) with A^ ' as in 

Eq.© and f) 2 0) = 1.667 x 1CT 3 . 
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di) irt) y/tf - (2/tt) 2 

Similarly to Eq. (flq) , this shows that as f) — > 2/tt, the two m 

solitons transform increasingly fast. 

Thus if we want to use the two real solitons (J43|) as ap- 
proximations for their respective (7 7^ 0)-counterparts, 
we should keep in mind that this approximation is valid 
only far away from the turning point 7thr- Since for 
small h the turning point is close to the /i-axis (i.e. 
7thr ~ (tt/2)/i is also small), the validity of the approxi- 
mation will be restricted to very small 7: j/h <C tt/2. 
2. Finally, we briefly summarize the main points of this 
work. 

(a.) The lower boundary of the existence domain of 
the two solitons is given by the following asymptotic ex- 
pression (as 7^0): 
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